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Abstract 


The reproductive behaviour of four intertidal and two subtidal mussel populations werc investigated within 
two spatial scales at Ninety Mile Beach, northern New Zealand, where large quantities of seed mussels are 
collected for aquaculture. Reproductive cycles were compared between intertidal and subtidal populations 
and among intertidal populations along a 90-km coastal area. Histological sections of mussel gonads were 
used to identify the reproductive state of individual mussels collected monthly between July 1998 and April 
2000. Gonad indices were calculated for populations cach month at each site. Frequency distributions of 
different gametogenic stages were plotted and differences in reproductive cycle among populations were 
identified. Mussel shell sizes were also recorded and compared among populations. Although all intertidal 
and subtidal populations had similar temporal variations in reproductive behaviour, spatial patterns among 
various sites revealed clear differences in reproductive cycles among populations. In particular, mussel shell 
size and gonad indices were higher for subtidal sites compared to their intertidal counterparts. Among 
intertidal mussel populations, the northernmost population had the highest maximum shell length, gonad 
indices and percentage ripe and spawning individuals, followed by two populations at the southern end of 
Ninety Mile Beach. A small mid-beach mussel population had the smallest mussels with the lowest gonad 
indices and few spawning mussels during the reproductive season. A prolonged spawning period was 
observed within the most reproductively active populations (two subtidal populations at the south end and 
one intertidal population at the north end of the beach) between June and December. However, the 
reproductive cycle of two intertidal populations at the south end of the beach indicated 2-3 spawning events 
during the same period. Differences in mussel shell size and reproductive behaviour among sites are likely 
to be related to variations in environmental conditions along Ninety Mile Beach. 


Additional keywords: gonad cycle, histology, intertidal, Perna canaliculus, spat, spawning synchrony, 
subtidal. 


Introduction 


The biology of the commercially important green-lipped mussel, Perna canaliculus, has 
received considerable interest in recent years because the industry largely depends on 
5160 tonnes per year of wild mussel spat (attached to macroalgae) collected from Ninety 
Mile Beach in northern New Zealand to seed their farms (Jeffs et al. 1999). However, 
understanding the basic biology of natural mussel populations is critical if mussel spat are 
to be utilised reliably in terms of their availability and arrival at coastal areas along Ninety 
Mile Beach. The reproductive behaviour of Perna canaliculus has been investigated for 
mussel-farm populations in Marlborough Sounds (Flaws 1975; Tortell 1976; Hickman et 
al. 1991; Buchanan 1998) and wild populations at Ninety Mile Beach (Alfaro et al. 2001). 
However, differences in environmental conditions, such as water temperature, salinity and 
food availability, make generalisations from one population to another difficult (Newell er 
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al. 1982: Hickman et al. 1991: Snodden and Roberts 1997). The temporal variability in 
reproduction of mussels at Ninety Mile Beach, where the majority of New Zealand's mussel 
spat is collected, has been investigated previously (Alfaro et a/. 2001). However, spatial 
variability in reproduction of mussel populations along the 90-km coastal stretch, as 
described herein, provides important additional information regarding the dynamics of 
gamete production in the area. These studies on the reproductive biology of Perna 
canaliculus are part of an ongoing, integrative research programme analysing larval 
behaviour, settlement and recruitment processes and adult mussel-bed population 
dynamics, which may contribute to the future sustainability of wild spat production from 
Ninety Mile Beach. 

Study of the reproductive behaviour of various mussel populations within a given 
geographical area can provide valuable information on variability of reproductive efforts at 
different spatial scales. This information then can be used to aid predictions of patterns of 
spatfall supplies to coastal areas, which can be utilised for commercial purposes and 
environmental resource management, One way to describe the reproductive efforts of a 
population is by using a gonad index as a standardised measure of reproductive maturity 
(Kautsky 1982; Brown 1988; King et al. 1989; Borcherding 1991: Galinou-Mitsoudi and 
Sinis 1994). Seasonal variations in reproductive cycle can be obtained and compared 
among populations by quantifying the frequency of mussels within various reproductive 
stages in the populations at a given time (Kautsky 1982; Shafec 1989; Villalba 1995; 
Nichols 1996). 

The aim of this study is to investigate differences in reproductive behaviour of Perna 
canaliculus populations within different habitats along Ninety Mile Beach, northern New 
Zealand. Reproductive patterns of four intertidal and two subtidal adult mussel populations 
from Scott Point to Ahipara Bay are examined and compared to address spatial variability. 


Materials and methods 


Six geographically distinct mussel populations likely to be affected by different environmental conditions 
were identified at Ninety Mile Beach, northern New Zealand. Four rocky, intertidal mussel populations 
were located at Tonatona Beach, Ungaunga Bay, The Bluff and Scott Point; whereas two subtidal 
populations were found at Wizard Rock and Blue House (Fig. 1). The Tonatona Beach and Ungaunga Bay 
populations are exposed to intermittent strong easterly swells from the prevailing easterly winds and 
northbound Westland current, as well as periodic sediment load disturbances from redistributed nearby 
sands. A mussel population at The Bluff (Fig. 1) is also exposed to both high wave action and sand load 
disturbances. The northernmost population at Scott Point experiences strong hydrodynamic regimes, but 
because of the extensive rock reef, the sediment load appears to be low in comparison with the other 
intertidal populations. The subtidal populations at Wizard Rock and Blue House are located ^8 km apart 
near Ahipara Bay and are in 13 and 15 m water depth respectively (Fig. 1). 

Gonad maturity analyses from histological sections were carried out on mussel specimens from the six 
study sites at Ninety Mile Beach. Approximately 70 adult mussels (50 to 145 mm in greatest length) were 
collected monthly from each of the six sites. An attempt was made to sample all sites once a month, but 
collections were not possible during certain months owing to poor weather conditions, or during months 
when collections were prohibited by cultural traditions of native Maori land owners, Sampling took place 
from July 1998 to April 2000. Intertidal mussels were collected haphazardly during low tide; whereas 
subtidal mussels were dredged when weather conditions permitted, from which 70 mussels were selected 
haphazardly for study from the dredge. Immediately after collection, the mussels were measured across 
their greatest shell length, sectioned dorsoventrally through the mid-body and the tissues fixed in Bouin’s 
solution. Mussel tissues were histologically processed and the resulting slides were analysed under a 
microscope according to techniques developed for this species (Alfaro ef al. 2001). Gonad indices, based 
on gonad maturity, for each population were calculated following King ef al. (1989) and Galinou-Mitsoudi 
and Sinis (1994) using the following formula: 
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(number in each stage) (rank of stage) - 


total sample size Së 


where the scoring system was: 


Resting stage 

Early development stage 
Late development stage 
Ripe stage 

Spawning stage 
Redevelopment stage 
Spent stage 


— NN LO b2 rs C5 


Frequency distribution graphs, including all developmental stages, were constructed to identify differences 
in gametogenic activity among the six sites investigated. 


Statistical analyses 


Mussel shell length and gonad index data were tested for any violations of the basic parametric statistical 
assumptions and no transformations were necessary. Repeated-measures designs, with a MANOVA 
approach and profile analysis methodology, were used on the mussel shell length and gonad index data to 
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Fig. 1. Map of the study site at Ninety Mile Beach, northern New 
Zealand. Four intertidal (Tonatona Beach, Ungaunga Bay, The Bluff 
and Scott Point) and two subtidal (Wizard Rock: and Blue House) 
mussel populations are shown. 
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test multivariate responses between intertidal and subtidal habitats and between southern and northern 
habitats. Thus, when comparing intertidal and subtidal habitats, two replicate intertidal sites and two 
replicate subtidal sites where used. Similarly, two replicate northern sites and two replicate southern sites 
were used in the southern v. northern habitat analyses. Using repeated-measures analyses made it possible 
to take irito account the correlation among times sampled (dates), which avoided the problem of 
pseudoreplication, because measurements among dates were not independent (i.e. the same population was 
sampled over time). For this type of analysis, habitat and date were always fixed factors and the sites nested 
within habitats were random. When the repeated-measures data were analysed with MANOVA, the 
response variable for each level of the within-subject factor (date) was treated as a different dependent 
variable, These analyses allowed us to test whether the mean response vectors among the between-subject 
(habitat) levels were different. In order to satisfy conditions of MANOVA, equal sample sizes were achieved 
by only using sampling events that were shared by all sites in a given statistical test. In addition, the large 
sample sizes of individual mussels ensured powerful statistical analyses. The pattern of response of the date 
factor was analysed with profile analysis, which uses the contrasts or differences of adjacent date-pairs to 
detect significance of treatments and interactions (Scheiner and Gurevitch 1993). Individual ANOVAs, 
with Bonferroni o adjustments, were also conducted for each significance test to identify specific dates 
when treatment effects differed. 


Results 
Mussel shell length 


Maximum shell length measurements of adult mussels from intertidal and subtidal sites at 
Ninety Mile Beach resulted in marked differences in mean maximum shell length among 
mussels from the six different sites. Generally, subtidal mussels were larger than intertidal 
mussels and mussels from Scott Point were consistently larger than mussels from all other 
intertidal sites. The Bluff had the smallest mussels of any site at Ninety Mile Beach 
throughout the study period (Fig. 2). Mean (+ s.e.) maximum shell lengths for Tonatona 
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Fig. 2. Mean (+ s.e.) maximum shell length of mussels within six adult mussel 
populations sampled monthly between July 1998 and April 2000. Means were 
calculated from 70 adult mussels on each date sampled within each site. Tonatona 
Beach, Ungaunga Bay, The Bluff and Scott Point represent four intertidal populations 
and Wizard Rock and Blue House constitute two subtidal mussel populations. 
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Beach, Ungaunga Bay, The Bluff, Scott Point, Wizard Rock and Blue House were 86 + 0.3, 
82 + 0.2, 68 + 0.2, 93 + 0.3, 101 + 0.4 and 99 + 0.4 mm, respectively; and ranges of shell 
length were 60-116, 59-111, 50-88, 64-125, 71-145 and 70-142 mm, respectively. 
Results from repeated-measures MANOVAs, comparing mean maximum shell length 
between intertidal (Tonatona Beach and Ungaunga Bay) and subtidal (Wizard Rock and 
Blue House) habitats and between southern (Tonatona Beach and Ungaunga Bay) and 
northern (Scott Point and The Bluff) habitats, are shown in Table 1. Single Bonferroni- 
adjusted paired contrast ANOVAs for variables, which revealed significant group 
differences, also are shown in Table 1. 


Female and male synchrony 


Perna canaliculus is dioecious and, from our results, no major differences were observed 
in gonad indices and reproductive cycles between female and male mussels in any of the 
six sites studied, nor in months sampled between July 1998 and April 2000 (Figs 4-5, 7). 
This high level of reproductive synchrony allowed for female and male gonad indices to be 
pooled for further comparisons among sites and dates. Temporal variability in reproductive 
behaviour was remarkably similar among all sites and is reported elsewhere (Alfaro et al. 
2001). Generally, populations showed higher levels of reproductive maturity and spawning 
between March and December and very little spawning outside this reproductive season. 
However, some differences in mussel reproductive behaviour were observed between 
habitats (Alfaro et al. 2001). 


Table 1. Statistical analyses of maximum mussel shell length 
A repeated-measures MANOVA (model III) was used to test multivariate responses between inter- 
tidal (Tonatona Beach and Ungaunga Bay) and subtidal (The Bluff, and Scott Point) habitats and, 
in a separate analysis, between southern (Tonatona Beach and Ungaunga Bay) and northern (The 
Bluff, and Scott Point) habitats. Profile analyses with contrasts between adjacent date pairs were 
conducted with individual ANOVAs and using Bonferroni o adjustments 


d.f. m.s. P p Significant 
> contrasts 
Intertidal v. subtidal habitats 
Between-subject 
Habitat 1 3628.15 52.42 0.019 
Error 2 69.22 
Within-subject 
Date 14 146.73 16.76 0.001 Jul/Aug, Aug/Sep, 
Habitat x date 14 19.60 eo Ba Sep/Oct, Oct/Dec 
Error (Date) 28 8.76 


Southern v. northern habitats 
Between-subject 


Habitat 1 1.57 0.001 0.973 
Error 2 1109.34 
Within-subject 
Date 7 17.57 Bal 0.032 
Habitat x Date 7 13.76 2.48 0.070 Jul/Aug, Aug/Sep 
Error (Date) 14 3:05 


Statistically significant contrasts between dates are shown on the right column (significant date 
pairs only occurred in 1999). 
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Intertidal v. subtidal populations 


The reproductive behaviour of two intertidal mussel populations (Tonatona Beach and 
Ungaunga Bay) was compared to the reproductive behaviour of two nearby (~10 km) 
subtidal populations (Wizard Rock and Blue House) at the southern end of Ninety Mile 
Beach. Intertidal mussels consistently had lower gonad indices than subtidal mussels 
(Fig. 3). The mean (+ s.e.) gonad index and range from Tonatona Beach mussels were 189 
+ 8 and 144-259, whereas Ungaunga Bay mussels had a gonad index mean (+ s.e.) and 
range of 173 + 6 and 136-247. The frequency distribution of all reproductive stages for 
mussels within Tonatona Beach and Ungaunga Bay indicated similar reproductive cycles, 
although at Ungaunga Bay, a greater number of mussels were in the spent stage throughout 
the post-spawning periods (Fig. 4). Mussels from both intertidal sites showed two or three 
likely spawning events between June and November that were concurrent at both sites, 
followed by rapid gametogenesis until the end of the spawning season in December (Figs 
3-4). Compared to the intertidal populations, the subtidal populations had gonad index 
means (+ s.e.) and ranges of 220 + 10 and 158—268 for Wizard Rock and 215 + 8 and 161— 
277 for Blue House, respectively. Furthermore, intertidal mussel populations had greater 
numbers of mussels in the spent stage of development than subtidal populations before the 
end of the spawning season in November/December (Figs 4-5). Conversely, the percentage 
of ripe and spawning mussels was consistently higher for subtidal mussel populations than 
intertidal populations throughout the reproductive season (March to December) (Figs 4-5). 

Repcated-measures MANOVA and profile analyses to test differences in gonad indices 
between intertidal (Tonatona Beach and Ungaunga Bay) and subtidal (Wizard Rock and 
Blue House) habitats among 15 months, between November 1998 and April 2000, are 
shown in Table 2. Repeated-measures ANOVAs to test differences in treatment effects 
between 14 paired-date contrasts also are shown in Table 2. 
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Fig.3. Gonad indices for two intertidal (Tonatona Beach, Ungaunga Bay) and 
two subtidal (Wizard Rock and Blue House) mussel populations at Ninety Mile 
Beach. A total of 70 mussels were sampled on each date and within cach site to 
calculate a gonad index at a given date and site. 
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Table2. Statistical analyses of mussel gonad index 
A repeated-measures MANOVA (model III) was used to test multivariate responses between intertidal 
(Tonatona Beach and Ungaunga Bay) and subtidal (The Bluff, and Scott Point) habitats and, in a sep- 
arate analysis, between southern (Tonatona Beach and Ungaunga Bay) and northern (The Bluff, and 
Scott Point) habitats. Profile analyses with contrasts between adjacent date pairs were conducted with 
individual ANOVAs and using Bonferroni a adjustments 


dt m.s. F E Significant contrasts 
Intertidal v. subtidal habitats 
Between-subject 
Habitat 1 19279.06 09132 0.042 
Error 2) 863.71 
Within-subject 
Date 14 3965.90 41.41 0.001 Oct/Dec, Dec/Feb, 
Habitat x Date 14 581.96 6.08 0.001 Feb/Mar, Apr/May, 
Error (Date) 28 95.78 Jun/Jul, Dec/Jan 
Southern v. northern habitats 
Between-subject 
Habitat 1 334.76 0.03 0.871 
Error 2 9881.36 
Within-subject 
Date 7 1694.58 5.17 0.004 
Habitat x Date 3j 357.98 1.09 0.419 Jul/Aug, Aug/Sep 
Error (Date) 14 327.76 


Statistically significant contrasts are shown on the right column (significant date pairs for intertidal 
v. subtidal habitats extend from 1998 to 2000, and occurred only in1999 for southern v. northern 
habitat comparisons). 


Southern v. northern populations 


Larger-scale comparisons in reproductive behaviour were conducted among southern 
(Tonatona Beach and Ungaunga Bay) and northern (The Bluff and Scott Point) intertidal 
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Fig. 6. Gonad indices for four interiidal mussel populations (Tonatona Beach, 
Ungaunga Bay, The Bluff and Scott Point) at Ninety Mile Beach. A total of 70 mussels 
were sampled on each date and within each site to calculate a gonad index at a given 
date and site. 
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mussel populations. Comparisons of gonad indices among the four intertidal populations 
indicated that The Bluff consistently had the lowest gonad indices (mean (+ s.e.) of 146 + 
4 and range of 136-177) during the sampled period, followed by the southern populations 
of Ungaunga Bay and Tonatona Beach (see previous section) (Fig. 6). Scott Point had the 
highest gonad indices of all intertidal populations, with a mean (+ s.e.) of 217 4 9 anda 
range of 180-262 (Fig. 6). Frequency distribution data suggest that the mussel population 
at Scott Point had the highest percentage of spawning mussels during the reproductive 
season (March to December) compared to southern intertidal populations, with The Bluff 
population showing the lowest percentage of spawning mussels (Fig. 7). Also, mussels in 
the spent developmental stage were more numerous in The Bluff population than the 
southern populations, whereas Scott Point had the lowest number of spent mussels during 
the reproductive season (Fig. 7). Repeated-measures MANOVA and profile analyses to test 
differences in gonad index between southern (Tonatona Beach and Ungaunga Bay) and 
northern (Scott Point and The Bluff) intertidal sites among 15 months between November 
1998 and April 2000 are shown in Table 2. Table 2 also shows the results of single ANOVA 
tests on paired-date contrasts. 


Discussion 


The timing and location of mussel spat arrival at Ninety Mile Beach have been previously 
unpredictable (Jeffs et a/. 1999). The complexity of interacting physical and biological 
parameters that affect the green-lipped mussel life cycle have made the study of the ecology 
of this mussel challenging in the past (Buchanan 1998). Research on the importance of 
various ecological processes of mussel populations at Ninety Mile Beach is underway and 
will likely provide evidence of the factors responsible for the high spatial and temporal 
variability of mussel spat production in the area (A. C. Alfaro, unpublished data). The study 
of spatial variability in reproductive behaviour among various intertidal and subtidal mussel 
populations contributes to a better understanding of the timing and magnitude of larval 
production at various locations along Ninety Mile Beach. 

Marked differences in maximum shell length among all intertidal and subtidal 
populations are likely a result of variations in environmental conditions that regulate 
feeding rates, such as food supplies, emersion time, hydrodynamic regime, sediment load 
and temperature (Suchanek 1985: Frechette and Grant 1991, Alvarado and Castilla 1996, 
Dahlhoff and Menge 1996, McQuaid and Lindsay 2000). Different environmental 
conditions may also result in different spawning times and duration, as well as more subtle 
variations in gametogenic cycles (Bayne ef al. 1983, Nichols 1996, Myrand et al. 2000). 
Experiments conducted by Myrand et al. (2000) in the Gulf of St Lawrence revealed 
distinctive differences in timing and duration of spawning events and in mortality between 
two populations of blue mussels, Mytilus edulis. The two mussel populations originated 
from the same initial stock but were subjected to different environmental conditions of 
temperature and water depth around the Magdalen Islands. In this study, comparisons 
between northern and southern intertidal populations and between intertidal and subtidal 
populations of Ninety Mile Beach resulted in differences in shell length between months 
throughout the year but site differences were only significant between intertidal and 
subtidal sites. Contrast analyses on single date pairs indicated that most of the shell length 
differences occurred between July and December 1999 for intertidal v. subtidal populations 
and between July and September 1999 for southern v. northern populations. However, the 
mussel shell length information provided in this study can only be used to recognise shell 
size differences among sites and.sampled dates. 
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The high level of reproductive synchrony observed between female and male mussels in 
this study has been reported in other Perna canaliculus populations (Flaws 1975; Tortell 
1976; Buchanan 1998) and in other mussel species (Heffernan and Walker 1989; Villalba 
1995; Nichols 1996; Darrigran et al. 1999). Despite the great physical distances among 
populations (from 8-90 km apart) at Ninety Mile Beach, a significant temporal synchrony 
was also observed in reproductive behaviours among populations over 90 km and both 
intertidal and subtidal populations (Figs 3-7). This synchrony may imply spawning 
triggered by some spatially reliable cue, such as changes in water temperature. 

These unique temporal patterns in mussel reproduction at Ninety Mile Beach indicate a 
prolonged spawning season between June and December, which differs from the two 
spawning periods reported for mussel populations in the South Island, New Zealand (Flaws 
1975; Buchanan 1998; Alfaro ef al. 2001). Conversely, the present study shows clear 
differences in gonad indices and reproductive cycles among mussel populations within 
different habitats at Ninety Mile Beach. 


Intertidal v. subtidal populations 


Comparisons of reproductive patterns among intertidal and subtidal mussel populations at 
Ninety Mile Beach indicate that subtidal populations have consistently higher gonad 
indices and a higher percentage of ripe and spawning mussels compared to intertidal 
populations. These trends in reproductive behaviour can be loosely associated with 
reproductive effort, although more detailed studies incorporating biomass and energy 
measurements are necessary. Subtidal populations have often been associated with higher 
productivity than intertidal populations (Suchanek 1978; Thompson 1979; Griffiths 1981; 
Reusch and Chapman 1997). Furthermore, subtidal mussels tend to have higher growth 
rates and reproductive efforts owing to the generally prolonged feeding periods and lower 
stresses from desiccation, temperature fluctuation and predation (Suchanek 1978; Griffiths 
1981; Bayne et al. 1983; Reusch and Chapman 1997). Thus, mussel populations at subtidal 
Wizard Rock and Blue House may sustain prolonged reproductive maturity by re-allocating 
energy from maintenance to reproduction. Another indication of higher reproductive 
activity is an increase in age and size, often associated with habitats that exhibit more 
constant environmental conditions (Griffiths 1981: Bayne et al. 1983: Petraitis 1995), such 
as those found at Wizard Rock and Blue House. Mussels from subtidal sites in this study 
not only had the highest synchronised gonad indices, but also maintained more ripe and 
spawning mussels through the season. Consistency of parameters in subtidal sites, such as 
temperature, compared to intertidal sites, also might account for the single but extended 
subtidal spawning period, with very few spent mussels in the populations until the end of 
the season. Conversely, rapid temperature changes in the southern intertidal populations at 
Ninety Mile Beach may be responsible for the distinctively rapid redeveloping and 
subsequent spawning periods through the season (March to December) for both of these 
populations (Alfaro et a/. 2001). The spatial and temporal reproductive differences between 
intertidal and subtidal populations may emphasise the importance of subtidal populations 
as major contributors to the spat production at Ninety Mile Beach. 


Southern v. northern populations 


Differences in mussel reproductive behaviour between intertidal habitats at Ninety Mile 
Beach suggest great variability in overall gonad indices among northern and southern 
populations. The mussel population at The Bluff had the lowest gonad indices. Although 
data for this site are not available from July 1999 to April 2000, there is an indication that 
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this population also had the shortest spawning period among all studied populations. It is 
difficult to explain why the population at The Bluff had such a low productivity when the 
nearest population (15 km to the north) at Scott Point had the highest reproductive output 
at Ninety Mile Beach. Indeed, Scott Point had a consistently high and long spawning period, 
with low numbers of spent mussels until November. It is possible that less favourable local 
conditions at The Bluff may contribute to the low reproductive productivity of this 
population. The Bluff is an isolated patch of rocks in the middle of a long stretch of sandy 
beach. Thus, the sediment load and wave exposure are high and may be a constant source 
of stress for mussels surviving on the rocks. Owing to the generally linear coastline of this 
area, food supplies tend to be irregular and may depend on the northbound nearshore 
current, with little residence time at The Bluff (A. C. Alfaro, unpublished data). Conversely, 
Scott Point is located at the northernmost end of the beach, where the mixing of the West 
Auckland and Westland currents, and the more complex coastal morphology of bays and 
shallow banks, together create a highly turbulent environment (Stanton 1973). Thus, the 
combination of dynamic ocean circulation and high delivery and residence time of food 
particles make the Scott Point area highly productive (A. C. Alfaro, unpublished data). The 
two southern populations appear to have intermediate reproductive behaviours compared to 
the Scott Point and The Bluff populations. Although statistical significance was only found 
for the date factor, contrasts tests indicated that differences between gonad indices were 
between July and September 1999. 

Regional variations in reproductive behaviour have been observed in many mussel 
populations throughout the world (Bayne et al. 1983; Thompson 1984; Jasim and Brand 
1989; Nichols 1996). This reproductive variability is believed to be a result of mussels 
responding to local conditions (Bayne et al. 1983; Nichols 1996). Thus, variations in 
environmental conditions such as water temperature and food availability may be sufficient 
enough to cause dramatic differences in reproductive behaviour among mussel populations 
(Hickman and Illingworth 1980; Smit et al. 1992; Dorgelo and Kraak 1993; Ram and 
Nichols 1993; Nichols 1996). For the green-lipped mussel, Perna canaliculus, great 
variations in condition index were found among populations throughout New Zealand, 
including populations ~36 km apart between Ahipara Bay and Whangaroa Harbour 
(Hickman and Illingworth 1980). Thus, the variations in reproductive behaviour among all 
intertidal and subtidal populations at Ninety Mile Beach may contribute, in various degrees, 
to the highly productive and sometimes spatially and temporally variable spatfall conditions 
observed in this unique coastal area. 

In summary, spatial variability in reproductive behaviour of mussels at Ninety Mile 
Beach were observed between different populations (between intertidal and subtidal 
populations and between northern and southern populations), possibly related to varying 
environmental conditions, such as temperature and food supply. Differences in reproductive 
behaviour between intertidal and subtidal populations may be a result of the more constant 
environmental conditions of subtidal populations. As a consequence, subtidal mussels tend 
to be larger (high growth rates) and have one prolonged spawning period. This spawning 
behaviour contrasts with the 2-3 spawning periods observed in the two intertidal 
populations and in Ahipara Bay. Regional differences in Ninety Mile Beach populations are 
exemplified by the disparate reproductive patterns among northern and southern mussel 
populations. Low reproductive maturity of mussels at The Bluff may be related to the high 
stress of sediment loads and irregular food supplies. Conversely, mussels at Scott Point 
benefit from higher water temperatures and high energy ocean conditions that may ensure 
high and constant food supplies to these mussel beds. These favourable mussel habitats may 
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result in the high reproductive activity of mussels at Scott Point and may promote a 
prolonged spawning season. Mussel populations at Ahipara have intermediate reproductive 
behaviours compared to those at The Bluff and Scott Point. 

Finally, the observed reproductive differences among mussel populations at Ninety Mile 
Beach suggest that the two subtidal populations at Wizard Rock and Blue House, and the 
northernmost intertidal population at Scott Point, had the highest reproductive output. This 
research examined thousands of mussels from an extensive 90-km reach of remote coastline 
in northern New Zealand. It is the first study to address the potential contribution of 
intertidal and subtidal mussel populations to the massive spat production of the area, which 
currently supports the commercial New Zealand mussel industry. Further studies on larval 
survival, transport and settlement will add important information regarding the temporal 
and spatial variability of spat arrival to the beach. 
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